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SIGNAL PROCESSING METHOD, SIGNAL PROCESSING SYSTEM, PROGRAM 
FOR SIGNAL PROCESSING, AND COMPUTER-READABLE STORAGE MEDIUM 
ON WHICH THIS PROGRAM IS RECORDED 

BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a signal processing method, 
a signal processing system and a program for signal 
processing, and a storage medium on which the program is 
recorded. In particular, the present invention relates to a 
signal processing method, a signal processing system, and a 
signal processing program that can be efficiently employed 
for a fast optical communication field, and a storage medium 
on which this program is recorded. 

Brief Description of the Prior Art 

Importance of fast and superior error correction technique 
In consonance with the expansion of the Internet and the 
development of e-business, the rate of increase in the 
volume of data computers can handle and their speed has 
accelerated. Accordingly, there is a demand for increasing 
speed of data transfer among computers, and in line with 
this demand, optical communication that yields transfer 
speed of up to 4 0 Gbps is becoming popular. However, for 
such a fast communication method such as the optical 
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communication method, to maintain an acceptable error rate 
at system level requires that the reliability of data 
communication be further increased in proportion to the 
amount of data processed by a computer. 

5 

Important techniques, called error correction coding 
techniques, have been devised to improve data reliability 
employing high-level mathematics to automatically correct 
errors caused by a variety of factors (e.g., noise along a 

10 channel) . Popular known techniques are Hamming code and 
Reed-Solomon code, both of which are frequently employed. 
Basically, Hamming codes correct single bit errors, but 
their correction capability is low. For instance, with 
Hamming codes, if a single bit error is detected, the error 

15 is corrected, but if the detected error covers two bits, 
only the error detection portion of the process is 
performed, no correction is made. However, setting up an 
error correction process for an error correction system 
using Hamming code is simple, and it is well known that by 

20 performing the error correction process in parallel, a 
processing speed that greatly exceeds 1 Gbps (one billion 
bits per second) can be obtained. 

On the other hand, Reed-Solomon code is a superior error 
25 correction technique possessing high correction capabilities 
and can be used to correct errors appearing as units 
(symbols) comprising multiple contiguous bits. At present, 
however, because of the complicated calculations that are 
frequently required, using Reed-Solomon code to perform 
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correction processes in parallel is difficult; and even when 
pipeline processes using data having eight-bit width are 
performed at 100 MHz, for example, only processing speed of 
around 800 Mb/s can be obtained. Currently, therefore, 
5 since the speed attainable with conventional techniques 
using Reed-Solomon codes is not suitable for fields in which 
high processing speed is required, these techniques are 
employed mainly for fields to which comparatively low data 
processing speed is acceptable, such as the low speed 
10 communication field and the data storage-unit field for the 
production of hard disks or secondary storage devices, 
CD-ROMs, for example. 

Error correction technique required by a fast optical 

15 communication field 

As part of a fast optical communication technique for data 
communication by computers, by the recently popular Internet 
as a backbone, a terabit per second fast communication 
system that employs WDM (Wavelength Division Multiplexing) 

20 and DWDM (Dense WDM), which has improved wavelength division 
multiplexing levels, has been introduced based on the SONET 
technique, according to which frames having a predetermined 
length are synchronously and sequentially transmitted. 

25 As the wavelength division multiplexing levels for the above 
described optical data communication is increased, crosstalk 
occurs between wavelengths that are near each other. To 
cope with this crosstalk, FEC (Forward Error Correction) is 
employed as an error correction method for communication in 
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long distance (Long Haul) optical wavelength division 
multiplexing. In ITU-T G.975, the ITU (International 

Telecommunication Union) standardized the use of interleaved 
(255, 239)RS code (code length n = 255 bytes) of m = 8 (8 
bits/symbols), and in G.709, the Digital Wrapper standard 
for defining the FEC frame structure is employed. 
According to the Digital Wrapper standard, for example, 
low-speed serial Reed-Solomon code circuits are arranged in 
parallel to achieve a necessary processing capability, and 
for this, interleaving of Reed-Solomon codes is an 
indispensable technique. 

Prior art for fast and high-level error correction 
techniques 

Independent of the need for optical communication, parallel 
fast decoding using Reed-Solomon code has been studied using 
a combinational circuit. 

Fig. 1 is a diagram showing an example of a fast decoder 
that can be used for an error correction apparatus. 
The decoding circuit in Fig. 1 implements a procedure for 
increasing by a multiple of ten or more the decoding speed 
of one decoding circuit, and for performing, at 
substantially the comparable speed as that of Humming code, 
parallel decoding in an error correction process using the 
Reed-Solomon coding possessing high-level error correction 
capabilities. For the decoding circuit in Fig. 1, a new 
representation using the elementary symmetric functions is 
employed for decoding Reed-Solomon codes, so that an error 
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value polynomial Er(x) of 0 (t) -degree with which error 
values can be directly calculated is employed (t is the 
maximum number of correctable errors) . 

Since the decoding circuit in Fig. 1 employs this 
polynomial, not only syndrome calculation and error location 
evaluation, but also error value evaluation can be directly 
obtained by performing calculations for a single polynomial. 

Therefore, compared with the conventional methods whereby 
these calculations are performed by Forney algorithm to 
divide evaluation results obtained for two polynomials, a 
greatly simplified operation can be used. Further, in the 
decoding circuit in Fig. 1, a representation appropriate for 
a combinational circuit is employed not only for the 
calculation of the coefficients of Er(x), but also for the 
calculation of the coefficients of the error locator 
polynomial A(x), so that not only can a higher processing 
speed be provided, but in addition, the number of required 
arithmetic circuits can be reduced. 

When the decoding circuit in Fig. 1 is employed, a random 
4-byte error correction circuit, which is mounted on an 
experimental base for a semiconductor using the standard 
0.35 pm ASIC technique, can process in parallel, and at a 
low latency (45 ns), data having a 320-bit width, and a 
processing speed of 7 Gb/s (7 billion bits per second) can 
be obtained that is nearly ten times higher than the typical 
processing speed of 800 Mb/s available with a current serial 
decoding circuit. In addition, when a new circuit 
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optimization algorithm specifically prepared for a large 
parallel error correction circuit and a circuit sharing 
method are employed for the decoding circuit in Fig. 1, the 
circuit size can be reduced. And furthermore, since the 
decoding circuit in Fig. 1 is a combinational circuit that 
does not require an external controller and registers, in 
spite of the high processing speed that can be attained, 
power consumption can be reduced. 

However, the decoding circuit in Fig. 1 can not provide a 
processing speed that equals 40 Gbps required for optical 
communication, and in order to cope with the 8-byte error 
correction standard established by the ITU, when the normal 
circuit sharing method is used the resulting circuit can be 
so large that it can not be mounted on a single chip. 

Fig. 2 is a schematic diagram showing the configuration of 
an error correction circuit that employs a conventional 
low-speed decoding method for optical communication. With 
this configuration, as the communication speed of an optical 
communication field increases, the conventional method 
whereby low-speed serial Reed-Solomon decoders are arranged 
in parallel becomes ever more inappropriate. Through 
conventional RS decoders have a processing speed below 1 
Gbps, the decoding method in Fig. 2 achieves the necessary 
processing speed by an appropriate arrangement of low-speed 
serial Reed-Solomon decoders. However, according to the 
conventional method in Fig. 2, for such an arrangement of 
many Reed-Solomon decoders are required, and accordingly, 
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the circuit size is increased in direct proportion to the 
data transfer speed used for optical communication. Fig. 3 
is a graph showing circuit size and data transfer speed 
plotted when the decoding method in Fig. 2 is employed. 

Fig. 4 is a diagram showing another conventional decoding 
circuit (A. Patel, IBM J. Res. Develop., vol. 30, pp. 
259-269, 1986) . Since according to this conventional 
decoding method, the processing speed can be easily 
increased for the calculation of syndromes and error 
locations. However, since as is shown in Fig. 4 Forney 
algorithm is employed for the calculation of error value, 
two polynomials, i.e., the differential dA(x)/dx for the 
error locator polynomial and the error evaluator polynomial 
Q(x), which are obtained by the syndromes and the error 
locator polynomial, must be evaluated, and then divisions 
must be performed. This is a critical path that prevents an 
increase in output speed, and the processing speed can not 
be satisfactorily increased . 

According to OC-768 SONET, this is a large problem, because 
assuming the 16 interleave defined by ITU-G709 is employed 
as an input/output interface for the decoding circuit, a 
fast processing speed of 300 MHz or higher is expected. 
Therefore, as one attempt, the decoder in Fig. 4 is employed 
and divisions corresponding to the critical path are 
converted into detailed pipelines to increase output speed. 

However, even when the process is converted into a pipeline, 



JP920010023US1 



the decoding circuit in Fig. 4 must perform divisions at 
locations whereat no error is present, and the circuit size 
and the power consumption are increased as the pipeline is 
constructed. Further, to perform divisions only for error 
locations, the locations must be calculated in advance, so 
that the error locations and the error values can not be 
calculated in parallel. In addition, for the decoding 
circuit in Fig. 4, a cycle count required for the output of 
the error values differs depending on whether an error is 
present. Therefore, when a synchronous frame, such as 
SONET, for sequential data must be input or output at high 
speed, it is difficult to output error values at high speed 
for a constant cycle, without depending on error patterns 
(number of errors and their locations) . 

Fig. 5 is a diagram showing an additional conventional 
decoder. When the parallel Reed-Solomon decoding method in 
Fig. 1 is employed for the optical communication field, 
because its circuit processing capability is superior to 
those of other conventional methods no problem occurs when 
non-interleaved RS code is used for an application. 
However, for interleaved Reed-Solomon codes, as defined by 
ITU-T G.975, since signals must be rearranged using a large, 
high-speed buffer and selector, the parallel Reed-Solomon 
decoding method is not always efficient. That is, the 
length of (255, 239) RS code is 2040 bits, and when a 16-byte 
interleaving process is performed, a 16-byte input and 
255-byte output serial/parallel converter and a 
parallel/serial converter for a 255-byte input and a 16-byte 
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output are required, thereby considerably increasing the 
size of a circuit even though the processing speed can be 
increased to a required level. Therefore, it is difficult 
for the parallel Reed-Solomon decoding method to be provided 
at a practical level for optical communication. 

For the calculation of error locations and error values used 
for the decoder, a large number of calculations in the 
Galois extension field GF(2 m ) must be performed at high 
speed, and further, the size of a circuit that can perform 
this processing must be such that it facilitates the 
implementation of the circuit. Conventionally, in the 
studies of the calculations over a Galois field, it is 
important that how efficiently a single calculation 
(multiplication or division) can be performed, and the 
several tens to hundreds of calculations by a combinational 
circuit have almost never been discussed to date. As one of 
various reasons this has not been done, it may be presumed 
many decoding operations tend to be performed by sequential 
circuits, and it has been ascertained that the use of a 
combinational circuit provides little merit in terms of 
processing capabilities and an acceptable circuit size. 

During the studies of the error correction calculation 
algorithm, the Yule-Walker equation that is defined for the 
Galois extension field GF(2 m ) is generated in decoding of the 
Reed-Solomon codes. The efficient processing of this 
Yule-Walker equation is desirable if high-speed processing 
is to be achieved and the size of the necessary circuit is 
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to be minimized. When the algorithm for solving the 
Yule-Walker equation is performed by a combinational circuit 
to achieve high-speed processing, in as the required error 
correction capabilities increase, the portion of the circuit 
used to solve the Yule-Walker equation and to locate errors 
becomes very significant from the viewpoint of the reduction 
in the size of the combinational circuit. 

In addition, when a combinational circuit that can carry out 
the decoding of Reed-Solomon codes is applied for an actual 
system, it is preferable that an algorithm be provided that 
can be applied for the decoding of Reed-Solomon codes having 
an arbitrary minimum distance in order to obtain a process 
that can be widely used and to remove superfluous, 
additional circuits or processes. Especially in the optical 
communication field, since the use of (255, 239) 
Reed-Solomon code is standardized by the ITU, an algorithm 
is required that can efficiently decode the Reed-Solomon 
code where the maximum number of correctable errors is 8 and 
the minimum distance is 17. 

In order to solve the mathematical problem posed by the 
Yule-Walker equation by using a hardware combinational 
circuit having a size that permits it to be implemented, 
increase in the circuit size must be suppressed, and an 
algorithm that can reduce the number of multipliers and a 
combinational circuit that can efficiently employ this 
algorithm are required. That is, a combinational circuit is 
needed that has an implementable size and that performs 
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high-speed processing, and that includes the error 
correction device and the error correction algorithm 
described above. 

SUMMARY OF THE INVENTION 

To resolve the above shortcomings, it is one object of the 
present invention to provide a signal processing method, a 
signal processing system and a signal processing program for 
efficiently processing interleaved Reed-Solomon codes in 
order to perform fast optical communication (40 Gbps or 
higher) , or more specifically, to perform SONET 
communication using wavelength division multiplexing, in 
which sequential data are transferred as synchronous frames, 
and to provide a storage medium on which the program is 
stored. 

That is, according to the invention, provided are a signal 
processing method and a signal processing system which have 
a high processing capability (low latency and high 
throughput) . 

It is another object of the present invention to provide a 
flexible signal processing method and a signal processing 
system that can process interleaved Reed-Solomon codes 
without the above described characteristics being lost. 

It is an additional object of the present invention to 
provide a decoding method and a signal processing system for 
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sequentially outputting error words at high speed and at a 
constant cycle rate, regardless of error patterns (error 
values and error locations) of the individual interleaved 
words that are received. 

It is another object of the present invention to provide a 
signal processing method and a signal processing system that 
can be used for a combinational arithmetic circuit, circuits 
over a Galois extension field GF(2 ra ) (m is an arbitrary 
natural number equal to or greater than 2), that has several 
common inputs and performs multiple multiplications (e.g., 
AB, AC and AD) and a circuit that performs logical sums such 
as AB + CD + EF + . . . can be implemented as small, fast, 
efficient processing circuits. 

It is a still further object of the present invention to 
provide a flexible program that can be applied for the 
decoding of Reed-Solomon codes having an arbitrary minimum 
distance and interleave configurations without an additional 
circuit or process being required, and to provide a storage 
medium therefor in order to perform the decoding of 
Reed-Solomon codes by a combinational circuit and to apply 
this combinational circuit to an actual system. 

The above objects can, according to the present invention, 
be achieved by providing a signal processing method, a 
signal processing system and a signal processing program, 
and a storage medium on which the program is stored. 
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More specifically, a signal processing method for a digital 
signal comprises the steps of: 

establishing a Yule-Walker equation having the 
following form by using a matrix that includes, as 
components, the elements of a Galois field GF(2 m ), and a 
vector that includes, as components, the elements of the 
Galois field GF(2 m ) 
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obtaining the solution of the Yule-Walker equation as 
the following determinants 





So 


Si 


" S1-1 




Ai = 


5 1-2-1 
Si-i+i 


Si-i ' 

Sl-i+2 ' 


S 21-i-2 

S21-1 


, i = i, . . . , 1-1 ; 




Si 


S1+1 


' S21-1 





employing Jacobi ' s formula, T ( / +1) Af t(i) +(Af t( ' } ) 2 =Ao at(i+1) ri , to 



enable the calculation of the solution (hereinafter 
referred to as Af t(i) ) to result in the calculation of the 
following determinants of the symmetric matrices 
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(where i - 0, . . . , 1) ; 

determining the number of errors to be the maximum 
matrix size that corresponds to the obtained solution that 
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is not zero; and 

determining whether the number of errors equals the 
maximum number of correctable errors. The components of the 
determinant are syndromes that consist of the elements of 
the Galois field GF(2 m ). The syndromes are generated by 
digital signals transmitted using wavelength division 
multiplexing. This invention is used for the decoding of 
digital signals or for error correction. 

According to the present invention, a system for processing 
a digital signal comprises: 

an encoding unit, for encoding a received digital 
signal ; 

a decoding unit, for decoding the digital signal that 
is encoded; and 

an output unit, for outputting the decoded digital 
signal, 

wherein the decoding unit includes 

means for establishing a Yule-Walker equation 
having the following form by using a matrix that includes, 
as components, the elements of a Galois field GF(2 m ), and a 
vector that includes, as components, the elements of the 
Galois field GF(2 m ) 
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means for obtaining the solution of the 
Yule-Walker equation as the following determinants 
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Ai = 
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Jacobi's formula, 
I^ 1+1) A^^+(Af*y=A^ I+1) ri r to enable the calculation of the 
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to result in 



solution Ai (hereinafter referred to as A* 
the calculation of the following determinants of the 
symmetric matrices 
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(where i = 0, . . . , 1) , 

means for determining the number of errors to be 
the maximum matrix size that corresponds to the obtained 
solution that is not zero, and 

means for determining whether the number of errors 
equals the maximum number of correctable errors. The 
encoding unit for encoding the received digital signals to 
syndromes that consist of the elements of the Galois field 
GF(2 m ). The syndromes are generated by digital signals 
transmitted using wavelength division multiplexing. This 
invention is used for the decoding of digital signals or for 
error correction. 



According to the present invention, a program for processing 
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a digital signal permits a computer to perform the steps of: 



'II 

Pi 



establishing a Yule-Walker equation having the 
following form by using a matrix that includes, as 
5 components, the elements of a Galois field GF(2 m ), and a 
vector that includes, as components, the elements of the 
Galois field GF(2 m ) 
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determining the number of errors to be the maximum 

matrix size that corresponds to the obtained solution that 
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is not zero; and 

determining whether the number of errors equals the 
maximum number of correctable errors. The components of the 
determinants are syndromes that consist of the elements of 
5 the Galois field GF(2 m ) . The syndromes are generated by 
digital signals transmitted using wavelength division 
multiplexing. The program of the invention is used for the 
decoding of digital signals or for error correction. 
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10 According to the present invention, a computer-readable 
storage medium is provided on which is recorded a program 
used for an error correction method, the program permitting 
a computer to perform the steps of: 

establishing a Yule-Walker equation having the 

15 following form by using a matrix that includes, as 
components, the elements of a Galois field GF(2 m ), and a 
vector that includes, as components, the elements of the 
Galois field GF(2 m ) 
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enable the calculation of the solution A i (hereinafter 
referred to as Af t(i) ) to result in the calculation of the 
following determinants of the symmetric matrices 

Si-i-i *•• S2(i-l-i) 52(1-2) "** ^21-1-i 
S l+l-i S 2 (l~i) ^2(2+1-2) 5 2 J+l-i 

Si "' S21-1-1 S21+1-X S21 
(where i = 0, . . . , 1) ; 

determining the number of errors to be the maximum 
matrix size that corresponds to the obtained solution that 
is not zero; and 

determining whether the number of errors equals the 
maximum number of correctable errors. The components of the 
determinants are syndromes that consist of the elements of 
the Galois field GF(2 m ) . The syndromes are generated by 
digital signals transmitted using wavelength division 
multiplexing. The program of the invention is used for the 
decoding of digital signals or for error correction. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a diagram showing a conventional decoder. 
Fig. 2 is a diagram showing a conventional error correction 
device for optical communication. 

Fig. 3 is a graph obtained by plotting a conventional 
circuit size and a data transfer speed. 

Fig. 4 is a diagram showing another conventional decoding 
circuit . 

Fig. 5 is a diagram showing an additional conventional 
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decoder . 

Fig. 6 is a schematic diagram showing a decoder according to 
one embodiment of the present invention. 

Fig. 7 is a diagram showing a multiplier having a 
conventional configuration . 

Fig. 8 is a diagram showing a multiplier having another 
conventional configuration. 

Fig. 9 is a diagram showing a multiplier having an 
additional conventional configuration . 

Fig. 10 is a diagram showing the embodiment wherein the 
present invention is applied for the multiplier in Fig. 7. 
Fig. 11 is a diagram showing the embodiment wherein the 
present invention is applied for the multiplier in Fig. 8. 
Fig. 12 is a diagram showing the embodiment wherein the 
present invention is applied for the multiplier in Fig. 9. 
Fig. 13 is a graph obtained by plotting the circuit size and 
the number of multiplier when the multiplier according to 
the invention is employed. 

Fig. 14 is a diagram showing a conventional error locator 
polynomial . 

Fig. 15 is a diagram of a formula, according to the 
invention, established for the Yule-Walker equation. 
Fig. 16 is a diagram showing the detailed structure of r ( 2 i+1) 
according to the invention. 

Fig. 17 is a diagram showing the detailed calculation 
results obtained, according to the invention, for decoding 
Reed-Solomon code . 

Fig. 18 is a schematic flowchart for an error correction 
algorithm according to the present invention. 
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Fig. 19 is a schematic diagram showing the configuration of 
a Reed-Solomon code decoder according to the invention. 
Fig. 20 is a schematic block diagram showing the 
configuration of an error correction device according to the 
invention. 

DETAILED DESCRIPTION OF THE INVENTION 
The preferred embodiment of the present invention will now 
be described by referring to the accompanying drawings. 
Note, however, that the present invention is not limited to 
this embodiment . 

Section 1 Decoding circuit 

Fig. 6 is a diagram showing a decoder according to the 
present invention that can be used to correct errors in 
digital signals received through optical communication. The 
decoder in Fig. 6 includes an input unit 10, a processor 12 
and an output unit 14. The input unit 10 receives a 16-byte 
interleaved digital signal ID. The processor 12 processes a 
signal received from the input unit 10 and calculates 
coefficients of an error locator polynomial and coefficients 
of an error value polynomial. And the output unit 14 
obtains an AND of the A {x) evaluation result and E r(x) 
evaluation result that are generated from the data received 
from the processor 12, performs an XOR process with the AND 
result and the input digital signal ID, and generates a 
digital signal OD. The digital signal OD is a signal in 
which errors that may have been included in the digital 
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signal have been corrected. 

In this invention, the digital signal ID input to the 
decoder in Fig. 6 can be one for which optical communication 
5 was used for its transmission, especially a signal 
transmitted by wavelength division multiplexing at a high 
data transfer rate of 40 Gbps. More specifically, the input 
digital signal ID can be, for example, a signal that is 
transmitted as (255, 239) RS code having a length of 2040 
10 bits. Generally, for wavelength division multiplexing, by 

jjsrls 

p using the interleave method, the input digital signal is 

y transmitted, for example, as a 16-byte parallel stream of 

pj 255 bytes to the decoding circuit of the invention. 

sas 

;; u 15 In the decoder design in the invention in Fig. 6, input 
digital signals ID are interleaved and are input in 
y parallel. A reception polynomial is defined for each input 

% signal, and a syndrome Si is calculated using this 

O polynomial and is output by the input unit 10. For (255, 

20 239) RS code, the syndrome Si output by the input unit 10 
(syndrome calculator) is generated as 16-byte digital 
information obtained from the 255 byte input digital signal. 
In the decoder for this embodiment in Fig. 6, the input 
digital signal of 2040 bits is interleaved, and the 16 
25 2 55-byte serial streams that are obtained are transmitted to 
the input unit 10, while sixteen 16-byte syndromes 
corresponding to the 16 serial streams are generated by the 
input unit 10 . 
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As is shown in Fig. 6, in the input unit 10, one syndrome 
calculator 16 is allocated for a serial stream IDS for each 
input digital signal ID, and calculates syndromes. The 
obtained syndromes for each serial stream IDS are stored in 
a register 18, and are output to the processor 12 one by 
one. For example, as is described above, 16 of 128 bit 
signals are obtained using the 16 interleaved 255 byte input 
digital signals in order to calculate error locations and 
error values. The syndrome calculator 16 that can be used 
for the input unit 10 in Fig. 6 can be any well known 
circuit, such as a circuit employing a sequential circuit. 
The method used for defining the syndromes and the 
calculation of the syndromes will be described later in 
detail . 

In the decoder in Fig. 6, the obtained syndromes Si are 
sequentially transmitted to the processor 12. For 
convenience sake, for the explanation only one processor 12 
is shown in Fig. 6. The processor 12 includes: an error 
locator polynomial calculator 18, which is a combinational 
circuit constituted by multiple multipliers and which is 
used to calculate an error locator polynomial A(x); and an 
error value polynomial calculator 20, which is used to 
calculate an error value polynomial Er(x). The 
combinational circuit that is used in the processor 12 and 
that includes multipliers will be described in more detail 
in Section 2 Combinational circuit>. 

The error locator polynomial A(x) and the error value 
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polynomial Er(x), which are output by the processor 12 in 
Fig. 6, are demultiplexed by a demultiplexer (not shown) so 
as to, for example, correspond to the number of interleaves, 
and the results are transmitted to the output unit 14. The 
output unit 14 includes registers 22, and AND gates 24 and 
XOR gates 26, which are arranged in a number equivalent to 
the interleaves of input digital signals. The output unit 
14 selects error values E r from the AND gates 24 by using the 
error location data A eva i ("1" for an error location or "0" 
otherwise) that are obtained by the syndromes. The selected 
error values E r are added by the XOR gates 26, and each XOR 
gate receives, via buffers 28a and 28b of imn bits, a serial 
stream IDS obtained from the coded digital signal ID. When 
the subtraction of the Galois extension field GF(2 m ) is 
performed by the XOR gate 26, the 255 byte digital signal OD 
from which an error has been removed is output. 

In the decoder in Fig. 6 for this invention, the 
combinational circuit constituting the processor 12 can be a 
sequential circuit . However, in this invention, especially, 
a plurality of multipliers may be constituted by three 
stages, an input side XOR calculator group (variable 
pre-processors) , an AND calculator group and an output side 
XOR calculator group (residual calculators) . When either 
one or both of the variable pre-processor and the residual 
calculator are used in common by multipliers, the size of 
the processor 12, which is the critical portion of the 
conventional calculation of error locations and error 
values, is practically allowable, and the processor 12 can 



JP920010023US1 



be efficiently constituted by multipliers. 

In the decoder of the invention in Fig. 6, the output unit 
14 is constituted only by circuits, such as a constant 
multiplier and an adder, performing linear calculations, 
without a circuit being used to perform non-linear 
calculations that reduce the processing speed. Therefore, 
for this invention a decoding circuit having a small circuit 
configuration can be provided that performs a process 
without deteriorating the processing speed. Further, 
through the careful study to provide this invention, the 
present inventors found that a decoding circuit that is more 
flexible, faster and smaller than a conventional circuit 
could be constituted by using a combinational circuit that 
includes multipliers having a specific configuration, and an 
algorithm that permits the combinational circuit to 
efficiently calculate error locations and number of errors. 

The method or the algorithm of the invention for the error 
locator polynomial calculation will, along with the 
combinational circuit, be described in detail later. An 
explanation will now be given for the function and the 
operation of the error value polynomial calculator 20 
included in the processor 12 of the decoding circuit of the 
invention . 

Selection of an error value calculation algorithm 

According to the invention, an algorithm that can directly 

evaluate not only error locations but also error values 
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through the calculation of the 0(t) -degree polynomial 
(linear calculation) is applied for the decoding of 
interleaved Reed-Solomon codes. At this time, in the 
algorithm used for this invention, the division necessary 
5 for an error value calculation is not performed for each 
error location on a critical path that is output following 
the evaluation of a polynomial, but by using only one 
calculation for each code word before the evaluation of a 
polynomial. Therefore, the values obtained by the 

10 polynomial evaluation can be directly output as error values 
Q at high speed at a constant cycle rate. 

v.ll 

W ST 

Hi Further, it has been found that the degree of the error 

p value polynomial Er(x) can be reduced to the t-1 degree that 

15 is the least required to obtain t independent outputs. At 
P this time, t coefficients can be calculated by using the 

coefficients of the error locator polynomial and the 
syndromes. By using this algorithm, not only the syndrome 
calculation and the error location evaluation, but also the 
20 error value evaluation can be performed merely by a linear 
operating circuit. Thus, the entire input/output circuit 
can be simplified and the processing speed can be increased. 

Various error value polynomials can be used for the decoding 
25 algorithm or the decoding method of the invention. In the 
explanation of the function and the operation of the 
invention, the following example is employed wherein when e 
errors have occurred at i 0 , . . . and i e -i the error value 
polynomial Er(x) is provided as 
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E 1JL E[(x+a J ) II (5 2 ^+a^) 

Er {e) (x)= 2 -i— r r r- 

v / 1=0 E[(a J +a *) 

and where, while the division is included, the denominator 
is not a polynomial but is a constant for each code word (in 
the above equation, a denotes the primitive element of the 
Galois extension field) . 



The error value polynomial Er (e) (x) can not be used as it is 
for the decoding circuit in this invention, because the k-th 
error value E ik (hereinafter referred to simply as E lk ) , 
p 10 which is present at location a lk (hereinafter referred to 



SI simply as a lk ) , is required for the calculation of the 

*P polynomial Er (e) (x) , and this makes no sense for the purpose 

O 

iy of this invention. However, if all error values E ik in 

* Er (e) (x) and locations a ik are written by the syndromes Si, 

h\ 15 this process can be provided as a circuit. Further, when 
all a ik s in Er {e) (x) are described by the error polynomial A^ e) 
, the error value can be calculated before the error 
locations are acquired . Therefore, the calculation of error 
values can be performed in parallel, and as a result, the 
20 processing speed can be increased. 



fl.j 



The above process can be performed by using the following 
decoding algorithm or decoding method. First, the 

denominator is written by A^ e) . The coefficients of the 
25 error value and error locator polynomials are elementary 
symmetric functions for error locations 

-jll o 2 2 
f CL f Q- f • • • i CL 

excluding a constant factor. For example, for coefficients 
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of the error locator polynomial, Ai, A 2 , . . . , A t , the 
following equations are established: 
Ax-Sa 1 * 

k 

A 2 = E a ik a^ 

k<l 

A t -. 1 = a i °a il a i2 . . .a 1 ^ 1 . 

The error locations a u r a il , a u , . . . , a 1 "- 1 can be exchanged with 
each other. Further, the following new relation 

m>l 

A (e) A (e) ... 0 ... 0 

1 A ( 2 e) Af -■ 0 -0 
0 A[ e) A ( 3 e) - 0 -0 

0 - 0 -• A { e % A ( e e) 
0 0 - A& A ( ^\ 

can be obtained by a relation (e.g., see "Symmetric 
Functions and Orthogonal Polynomials", I . G. Macdonald, 
American Mathematical Society, 1998) that is established 
between the Elementary symmetric function and the Shur 
function defined by the division of two Vandermonde 
matrixes, and the exchange relation, (a + b) 2 = (a - b) 2 = a 2 
+ b 2 , that is established between the addition and the square 
performed for the Galois extension field GF(2 m ). 

When the above determinant is employed, the denominator of 
Er (9) (x) can be written by A ( 2 e) . When the numerator of Er {e) (x) 
is calculated in the same manner, as is indicated by the 
following equation, the coefficients of Er (e) (x) can be 
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written merely by Si and Ay , and can be calculated without 
using E lk and a ik . 

Er^ e \x) 

lE ik U(x+a^) II (a^ + a 1 *) 

S 1 E- r£^A (e) }f(*-V(A {e) A (e) A (e) ^ 

= f^A^A^, ...,A ( e e) ) 

Z B lt a^Brg(A« f A ( 2 e) , . . . , A<?>) 

f^A^A^, ...,A ( e 6) ) 

.1 S m xlEr^(Kf,kf f ...,A ( e e) ) 
fW(A«,A» f ...,A ( e e) ) 

where A^^A^+a^A^I^ . Specifically, this coefficient, is 
that of the error locator polynomial that corresponds to the 
error location except a lk . 

As is described above, when these new relations are 
employed, the output unit 14 can be constituted by a linear 
operating circuit for the Galois extension field GF(2 m ), 
without a non-linear operating circuit being required, and 
the processing speed can be increased. 

The following configuration is employed to constitute a 
fast, small coding circuit that can also cope with 
interleaved codes based on the above error value calculation 
algorithm. 
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(1) Use of fast input unit and output unit (linear operating 
circuit) 

First, in consideration of the structure of the code (the 
number of interleaves) , the bus width of an input/output 
5 interface and the number of clocks, a polynomial evaluation 
(linear operating) circuit, which is connected to the input 
side to calculate syndromes, and a polynomial evaluation 
(linear operating) circuit, which is connected to the output 
side to evaluate error locations and error values, are 
M 10 implemented as fast sequential circuits that employ a cyclic 
structure that comes from the fact that RS codes are cyclic 



O 



JO codes, and that perform a pre-process and a post-process for 

r (n, k) Reed-Solomon code at an arbitrary number of clocks 

from 1 to n. Especially effective for the invention is the 
5 15 fact that the above described structure is also employed for 
an error value evaluation portion that is a critical path 
Q for the conventional method. According to the structure of 

% the invention, a decoding circuit can be provided as an 

pj interface that can flexibly cope not only with an input 

20 digital signal that has 255 bytes for each code word, but 
also with ones that have 1, 3, 5, 15, 17, 51 and 85 bytes. 

Table 1 shows the relation between the width of a digital 
signal input to the decoding circuit and the processing 
25 clock required for this input width. 

Table 1: The input/output width and a clock count for each 
code word (n = 255) 
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Input /output 
width (bytes) 


1 


3 


5 


15 


17 


51 


85 


255 


Clock count 
required for 
processing a code 
word (n = 255) 


255 


85 


51 


17 


15 


5 


3 


1 



As is shown in Table 1, although the required number of 
clocks is increased as the width of the input digital signal 
is reduced, the decoder of the invention can flexibly cope 
with it. 

In addition, in this invention, an arbitrary input/output 
byte width other than those in Table 1 can be selected for 
each code word. For example, an input/output width of eight 
bytes can be coped with when the code length is n = 256 
bytes by adding a one byte dummy at the end. 

(2) Connection to a non-linear operating circuit 12 
In this invention, a plurality of input units 10 and output 
units 14 that are constituted by sequential circuits are 
prepared in a number equivalent to the interleaves. Between 
these units, the non-linear operating circuit 12 is 
connected that calculates the coefficients of an error 
locator and error value polynomials by a multiplexer, a 
demultiplexer and a circuit for holding the syndromes and 
the coefficients of an error locator and error value 
polynomials. As is described above, this operating circuit 
12 is constituted as a combinational circuit for a 
non-linear operating circuit, such as a multiplier that 
performs non-linear calculation. Therefore, in this 

invention, in the case where ITU-T G.709 (255, 239) 
Reed-Solomon code is obtained by 16-byte interleaving, 16 
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sequential circuits must be prepared before and after the 
operating circuit 12. In other words, in this case, 
multiplexing and demultiplexing must be performed at a ratio 
of 16:1. However, in this invention, in order to multiplex 
5 and demultiplex the coefficients of the syndrome polynomial, 
only a signal of 128 bits (or 136 bits at the succeeding 
state, depending on how A(x) is defined) need be processed 
instead of 2040 bits, so that the required numbers of 
multiplexers, demultiplexers and buffers can be greatly 

1^ 10 reduced. 

Q 

fee* 

ji (3) Three-stage pipeline operation system that employs the 

non- linear operating circuit in the center in a time sharing 
fi manner 



s 

Q 



U 

;; ■: ; 

n 



15 According to the present invention, the entire decoder is 
operated as a three-stage pipeline formed from a linear 
operating circuit (syndrome calculation), a non-linear 
operating circuit (calculation of coefficients for error 
locator and error value polynomials) and a linear operating 

20 circuit (evaluation of error locations and error values), 
and the non-linear operating circuit of a low latency is 
employed in a time sharing manner by serially providing 
syndromes for code words for the calculation of interleaved 
code. Therefore, it is possible to provide efficient 

25 decoding of interleaved codes whose processing capability 
per circuit size is high. In the OC-768 case, for example, 
in order to operate the three-stage pipeline, the non-linear 
operating circuit in the center must complete the process 
for each interleaved code word at a latency of about 40 ns. 
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When the combinational circuit is implemented by the most 
advanced semiconductor technique (0.18 \im or better), the 
above decoder can be provided as an ASIC semiconductor 
device . 

More specifically, in the invention, fast decoder and error 
correction devices, the sizes of which fall within an 
acceptable range, can be provided by synergistic effects 
obtained by: (1) the use of a fast input unit 10 and a fast 
output unit 14, especially the use of only a linear 
operating circuit when calculating error locations and error 
values , the application of the above described algorithm for 
the processor 12, and the use of a configuration for a 
multiplier that can reduce the circuit size, and (2) the 
employment of the three-stage pipeline operation system by, 
in a time sharing manner, the processor 12, which 
essentially performs a non-linear calculation. In section 
2, which follows, an explanation will be given for a 
combinational circuit that comprises multipliers included in 
the processor 12 of the decoder of the invention. 

Section 2 Combinational circuit 

The processor 12 used for the decoder of the invention is a 
nonlinear circuit, specifically, a combinational circuit 
that uses multipliers. Unlike multipliers used for a 
conventional combinational circuit that has two stages, 
consisting of an AND calculator group and an XOR calculator 
group, in order to perform the multiplication of the Galois 
extension field GF(2 m ), the multiplier used for this 
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invention has three stages, consisting of XOR gates, AND 
gates and XOR gates. 

Configuration of a single parallel-multiplier 

While many studies have been made of a single multiplier, a 
parallel multiplier (Mastrovito Multiplier) , which is 
constituted as a combinational circuit, not as a sequential 
circuit, is a field of active research. For a conventional 
parallel multiplication circuit (hereinafter, in this 
specification, referred to simply as a multiplier) , there 
are two configuration types: the AND-XOR type and the 
XOR-AND-XOR type, which can be converted into each other. 
It should be noted, however, that the AND-XOR type is 
generally employed when a circuit is provided for only a 
single multiplication. This is because while the AND-XOR 
type has been well studied and various methods for obtaining 
a small circuit have been proposed, there is no guarantee 
that the circuit size for the XOR-AND-XOR type will be 
reduced (or may be increased) , and the reduction effects can 
not be obtained that would compensate for the expenditure of 
the effort a complicated design operation would entail. The 
AND-XOR type and the XOR-AND-XOR type will be further 
described. 

(1) AND-XOR type 

This is a typical method used for performing calculations in 
the same manner as are calculations performed with figures 
written down on paper, and generally a circuit of this type 
is employed. Specifically, the coefficients of two (m-l)th 
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degree polynomials, which are multiplication arguments, are 
combined to prepare m 2 partial products. This is the 
processing performed by the AND section. Then, partial 
products that have the same degree are added together to 
form a (2m-2)th degree polynomial, and a residue operation 
using an irreducible polynomial is performed to obtain the 
(m-l)th degree solution. This is the processing performed 
by the XOR section. The number of ANDs is m 2 and the number 
of XORs is 0(m 3 ), and it is widely known that (m 2 - 1) XORs 
can be obtained by selecting an appropriate irreducible 
polynomial and basis. An arbitrary multiplier can always be 
constituted by this method. 

(2) XOR-AND-XOR type 

Generally, according to the Boolean algebra rule (A and B) 
xor (A and C) = A and (B xor C) , it is possible for the XOR 
calculation in the residue operation unit of the AND-XOR 
circuit to be moved in front of the AND for use as a 
variable pre-processor (input side XOR calculator) . Thus, 
the XOR-AND-XOR multiplication circuit can be obtained. For 
moving the XOR calculation, when an even number of the same 
redundant terms are added to the XOR of the residue 
operation unit (output side XOR calculator) by the 
properties A xor A = 0 and B xor 0 = B, many XOR 
calculations may be moved to the front-end as a 
pre-processor . Since this operation can be employed, 
various methods for moving the XOR are available, in 
addition to the simple application of the distributive law. 
Therefore, multiple AND-XOR types are present even with the 
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same basis or irreducible polynomial. The number of gates 
in the XOR-AND-XOR type can vary, and may become either 
greater or smaller than that in the AND-XOR type. Another 
method is known whereby the number of XOR gates is 
systematically reduced by choosing a special basis, such as 
the Composite Field Multiplier that will be described below. 

(3) A Method for constructing an XQR-AND-XOR type that can 
be applied only for a limited field (Composite Field 
Multiplier) 

The composite field multiplier is a multiplier construction 
method that can be used only in a special case, such as 
where m is a composite number and the basis used for the 
representation of an element in a field may not be an 
ordinary basis (such as polynomial basis or normal basis) . 
This method will now be described in detail. When m is a 
composite number, extension field GF(2 m ) can be constructed 
by two fold extensions of the field GF(2). The composite 
field multiplier method is a method for constructing, in 
accordance with the extension process, a multiplier having a 
recursive structure. At this time, when the product of the 
two values Ax + B and CX + D in GF(2 m ) (A, B, C and D are 
values in the sub-field GF(2 m/2 )) 

(Ax + B) (Cx + D) = ACx 2 + (((A+B) (C+D) + AC + BD) x + BD) 
is employed for one quadratic extension, the number of 
multiplications performed in the sub-field can be reduced 
from four to three, and the circuit size can be reduced 

(KOA) . At the same time, the circuit can be provided as an 
XOR-AND-XOR structure (the addition performed before the 
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multiplication corresponds to the XORs arranged in front of 
AND) . It should be noted that the use of KOA is premised on 
the use of the composite field multiplier, otherwise, KOA 
can not be used. If the value of m is the composite number 
for which this method can be applied, a field converter is 
required, and the circuit size is increased because of the 
overhead. Thus, when a circuit is to be prepared for only a 
single multiplication operation, the circuit structure does 
not generally, employ the composite field multiplier. 

Conf iguration of a combinational circuit according to the 
invention using ordinary multipliers 

Figs. 7 and 8 are diagrams showing a configuration wherein 
multipliers with common input and a logical sum calculation 
circuit for the invention are constituted by a common 
AND-XOR type. In the example in Fig. 7, a combinational 
circuit using two multipliers is shown. As is shown in Fig. 
7, a first input Al is transmitted to both a multiplier 40 
and a multiplier 42, while a second input Bl is transmitted 
to the multiplier 40 to obtain a first output 45, and a 
third input B2 is transmitted to the multiplier 42 to obtain 
a second output 46. In the conventional structure in Fig. 
l f even though the multipliers receive common input, no 
circuit exists that these multipliers can employ in common. 
In the example in Fig. 8, the conventional multiplier 
structure is employed for a combinational circuit that 
performs logical sum calculations. It is again apparent 
that in Fig. 8 no circuit exists that the multipliers can 
use in common when performing logical sum calculations. 
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Fig. 9 is a diagram showing a combinational circuit 
employing conventional multipliers. In Fig. 9, one symbol 
is represented by one line, and 8-bit width input and output 
are assumed. The combinational circuit in Fig. 9 includes 
six-symbol input and one-symbol output, and seven 
multiplication circuits and five addition circuits. In the 
combinational circuit in Fig. 9, inputs SO, . . . and S3Q 
are transmitted to multipliers 46 and are added together by 
adders 48. The results are then transmitted to a logical 
sum calculation circuit 50 whereat an output L21Q, the 
logical sum of the inputs, is generated. 

In Fig. 9, a combination of cross-term generation operation 
and residue operation, indicated by broken lines, 
corresponds to one multiplication. Since a standard circuit 
structure is used for the multipliers shown in Fig. 9, no 
detailed explanation will be given for the circuit. The 
combinational circuit, including the multiplication circuit, 
includes 64 AND gates and about 103 XOR gates, and for the 
entire circuit, the number of gates is 448 ANDs + about 761 
XORs. As is apparent from Fig. 9, for almost all the 
multipliers, one or both inputs are used in common by 
another multiplier. Further, the logical sum calculation is 
performed at the final stage. 

Table 2 shows the numbers of gates that are included in 
multiplication circuits for a Galois extension field GF(2 8 ). 
They are conventional multiplier which has two-stage AND-XOR 
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structure, a composite field multiplier, and a 
multiplication circuit, for which XOR-AND-XOR is employed 
for the multiplication of a sub-extension field GF(2 4 ) . 



5 Table 2 



M 

J; 
13 



01 

o 



Variable 
pre-process 

None 



Standard AND-XOR 
multiplication circuit 

Composite Field Multiplier 4XOR*2 



Cross-term 
generation 

6 4 AND 
4 8 AND 



Residue 
operation 

+ 103XOR 
+ 56XOR 



XOR-AND-XOR circuit for performing the multiplication of sub-extension field 
GF(2 4 ) to obtain the effects of the invention 



a. Alteration of all three 22XOR*2 
multiplications of a 

sub-field 

b. Alteration of two out 16XOR*2 
of three multiplications 

c. Alteration of one of 10XOR*2 
three multiplications 



30 AND 



36AND 



42AND 



4 4XOR 



4 8XOR 



52XOR 



10 As can be seen from Table 2, as far as only a single 
multiplier is concerned, the circuit sizes for cases a. to 
c ,. are increased when compared with the size of the 
conventional composite field multiplier, and are not the 
minimum size. As is described above, when the three-stage 

15 XOR-AND-XOR structure is simply employed for a single 
multiplier, in some cases the circuit size may be increased. 



Multiplier structure for a combinational circuit according 
to the invention 
20 Generally, optimization of Boolean algebra is difficult when 
many multiplications and logical sum calculations are 
performed together. However, in consideration of the use of 
the combinational circuit as the processor 12 for the 
invention, many logical sum calculators and many multipliers 
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are connected in parallel at multiple stages. The i-th 
logical sum calculator generally receives 0-th to (i-l)th 
outputs. Therefore, since the backend part of the operation 
circuit must process in parallel inputs used in common by 
almost the entire combinational circuit, the optimization 
range is extended. The present inventors focused on this 
point, and achieved an efficient configuration for a 
multiplier by performing the optimization for Boolean 
algebra, while taking into account the balance obtained with 
other operations. 

Fig. 10 is a diagram showing a combinational circuit 
according to the embodiment wherein a three-stage structure 
is employed for the multiplier in the combinational circuit 
in Fig. 7, which includes conventional multipliers and 
adders. In the combinational circuit in Fig. 10, an input 
Al is transmitted to a first XOR group 52, an input Bl is 
transmitted to a second XOR group 54, and an input B2 is 
transmitted to a third XOR group 56. These XOR groups 52, 
54 and 56 are gates used by this invention to perform the 
variable pre-process. Since the XOR groups 52 process the 
common input Al, the circuit size is reduced. The outputs 
of the XOR groups 52, 54 and 56 are transmitted to AND 
groups 58, and when a cross-term is obtained, the residue 
operation is performed again by downstream XOR groups 60, an 
output 62 being generated by an XOR group 60a, and an output 
64 being generated by an XOR group 60b. In Fig. 10, the 
unit for performing one multiplication is indicated by 
broken line BL, and three stages, a variable pre-processor 
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(XOR) - a cross-term operation unit (AND) - a residue 
operation (XOR) , constitute one multiplier. As is shown in 
Fig. 10, when the XOR-AND-XOR structure is employed for each 
multiplier, and when for the multipliers the XOR 
calculations to be performed for an input used in common is 
unified, the XOR calculation circuits can be shared by the 
multipliers. So long as the size of a portion (variable 
pre-processors*2 + new cross-term generation + new residue 
operation unit) that corresponds to one multiplier is the 
same or slightly larger than the normal multiplication 
circuit, the size of the entire multiplier group can be 
reduced. 

Fig. 11 is a diagram showing a combinational circuit 
according to another embodiment wherein a three-stage 
multiplier used for the invention is mounted for the 
conventional combinational circuit in Fig. 10. In the 
combinational circuit in Fig. 11, inputs 66, 68, 70 and 72 
are transmitted to XOR groups 74, 76, 77 and 78 that perform 
the variable pre-process, and the outputs of the individual 
XOR groups are transmitted to AND groups 80 and 82 that 
perform the cross-term generation operation. 

The outputs of the AND groups 80 and 82 are transmitted to 
an addition circuit 84 and are added together, and the 
residue operation and the multiplication are again performed 
by a backend XOR group 86 that is used in common, so that an 
output 88 is generated. The structure of one multiplier in 
Fig. 11 is formed in a block Bx. The difference from the 
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structure in Fig. 10 is that, in the invention, even when 
the XOR groups 74, 76, 77 and 78 on the input side are not 
shared the backend or output side XOR group 86 that performs 
the residue operation is used in common. 

Further, in the invention, when the XOR groups 74 to 78 on 
the input side, which perform the variable pre-process, are 
used in common, and the XOR groups that perform the residue 
operation are also used in common, the overall structure of 
the combinational circuit can be simplified. 

Fig. 12 is a diagram showing a combinational circuit 
according to an additional embodiment wherein a three-stage 
structure is employed for the multiplier in the 
combinational circuit in Fig. 9. The combinational circuit 
in Fig. 12 constitutes a part of the processor 12 in the 
decoder (e = 2) in Fig. 6 for error correction by RS codes. 
The combinational circuit in Fig. 12, as does the 
conventional example in Fig. 9, receives input SO to S3Q, 
and as is shown in the blocks in Fig. 12, an XOR gate 90 
that performs the variable pre-process for input SO is used 
in common by three multipliers that correspond to AND gates 
92, 94 and 96. A residue operation unit 98 is shared by a 
plurality of multipliers, and both an XOR gate used for the 
variable pre-process and an XOR gate used for the residue 
operation are used in common. One multiplier is indicated 
by a chained line. In the combinational circuit in Fig. 12, 
eight variable pre-processors , seven cross-term generation 
operation units, four residue operation units, two 8-bit 
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adders, and three adders having the same bit width as the 
cross-term generation operation unit are provided. 

Therefore, when, based on the multiplication circuits in 
Table 2, a part of the decoder is to be constituted by the 
combinational circuit in Fig. 12, the number of gates 
required can be determined as shown in Table 3. 



Table 3 

Composite Field Multiplier 
Present invention 

a. Alteration of multiplication for all 
three sub-fields 

b. Alteration of multiplication for two of 
the three sub-fields 

c. Alteration of multiplication for one of 
the three sub-fields 



416XOR + 336AND Total 752gate 

458XOR + 210AND Total 668gate 

444XOR + 252AND Total 696gate 

430XOR + 294AND Total 724gate 



It is apparent from a. to c. in Table 3, that even if the 
single multiplier is not intentionally minimized, the 
overall size of the circuit can be reduced. 



Fig. 13 is a graph showing the number of XOR gates and the 
number of corresponding multipliers that are required for 
use for the processor 12 in the error correction circuit 
having the error correction capability t = 2 to 8 in Fig. 6. 

In Fig. 13, the vertical axis represents the total number 
of XORs, and the horizontal axis represents the number of 
multipliers. In this case, for the decoder for error 
correction in Fig. 6, m = 8 and the irreducible polynomial 
is x 8 + x 4 + x 3 + x 2 + 1. In this case, 662 multiplications, 
531 additions and 30 square calculations are performed. As 
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is apparent from Table 3 and Fig. 13, even if the single 
multiplier is not intentionally minimized, by the 
configuration of the invention, a greater reduction can be 
obtained in the overall size of the combinational circuit. 

In the embodiment in Fig. 13, the same variable 
pre-processor , cross-term generation operation unit and 
residue operation unit were used for all the multiplications 
in order to clearly present a specific circuit structure. 
In order to obtain better results when these circuits are 
actually implemented, the number of XORs in the arrangements 
for the variable pre-processor and the residue operation 
unit (i.e., which of a. to c. is to be employed) can be 
changed and optimized for each multiplication operation 
performed by the circuit. Further, in the embodiment in 
Fig. 13, since the ratio of the number of operations to the 
number of inputs is small and an overhead for field 
conversion is also present, no very great effect is obtained 
by reducing the number of gates. However, in actuality, 
since as is shown in Fig. 14 the ratio of the operations to 
the input is quite high, because of the effects produced by 
a reduction in the number of gates the circuit size can be 
dramatically reduced. 

Section 3 Error correction algorithm 

A detailed explanation will now be given for an error 
correction algorithm that is used by the decoding circuit 
and the error correction device of the invention. 
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Overview of the conventional example 

A. Conventional method for solving the Yule-Walker equation 
or for obtaining an error locator polynomial, and a problem 
associated with these methods 
5 According to the invention, it is necessary to find an 
efficient algorithm for a combinational circuit to calculate 
the following simultaneous linear equation, which is defined 
over GF(2 m ) , 



h 



. & 

flj 

Q 



10 



So 
Si 



Si 

s 2 



\ Si_i Si 



Si-i 
Si 



S21-2 J 



A (i > J 



\ -0 21-1 J 



where S 0 , Si, . . . S 2t -i are the elements of a given GF(2 n 
and A { ? is an unknown amount. 



ill 



In this simultaneous linear equation, the matrix on the left 
15 has a regular structure wherein the same elements are 
arranged obliquely to the right (the direction intersecting 
the diagonal line), and is called a Hankel matrix. 
Generally, this type of equation is called a Yule-Walker 
equation, and it is known that this equation is widely 
20 applied for various fields, such as the error correction 
code theory, time-series analysis and signal processing. In 
the error correction algorithm, the Yule-Walker equation 
appears in the portion for determining an error locator 
polynomial. Therefore, according to the present invention, 
25 the algorithm for obtaining the solution for the Yule-Walker 
equation is applied as an error correction algorithm for 
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decoding the Reed-Solomon codes. 



The well known methods for solving the Yule-Walker equation 
are, for example, the algorithm proposed by Levinson and the 
5 algorithm proposed by Levinson-Durbin, These algorithms 
start calculating with a matrix of the smallest size 1 and 
recursively determine the solution of an equation wherein 
the size of a matrix is greater. The number of calculations 
required by these two algorithms is of order of l 2 . However, 
10 these algorithms include divisions in the calculation step. 
y k This means that when the algorithm is mounted as a 

combinational circuit, conditional branches occur depending 
on whether the denominator is 0 or not. Since a separate 
circuit must be prepared for each conditional branch, the 
15 required circuit size increases as the size of the matrix 
increases . 



ft ;, 
4-- 



fl 

W 

en 



Further, as the object of the invention, especially relative 
to the decoding of the Reed-Solomon codes, an error locator 

20 polynomial is determined by obtaining the solution to the 
Yule-Walker equation. The conventional methods for solving 
the Yule-Walker equation can be, for example, the Peterson 
method, the Berlekamp-Massey method and the Euclid method. 
These methods are used to calculate the coefficients of an 

25 error locator polynomial by calculations of which the number 
is polynomial order with respect to the maximum number of 
correctable errors t. However, when the Berlekamp-Massey 
method and the Euclid method are represented by a 
combinational circuit, the following problems occur. 
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First, for the Berlekamp-Massey method, it is inevitable 
that multiple conditional branches should be included in the 
algorithm. Therefore, to expand this algorithm for a 
5 combinational circuit, for the same reason as described 
above, the circuit size would be increased in accordance 
with the number of combinations. As for the Euclid method, 
while the multiplication and division of the polynomial are 
the essential part of the algorithm and, the degree of a 
10 polynomial that appears in the denominator of the division 
^ can not be identified in advance, so that there is room for 

fear 

p generating a conditional branch. Furthermore, due to the 

^ conditional branches, the circuit size is accordingly 

,JE increased, as it is for the Berlekamp-Massey method. 

H 15 

3 B. Policy for calculating the Yule-Walker equation and an 

m error locator polynomial that is appropriate for a 

O combinational circuit 

Since, as is described above, the Levinson (-Durbin) method, 

W 20 the Berlekamp-Massey method and the Euclid method include 
conditional branches, a problem has arisen in how to provide 
these methods as combinational circuits. In order to 
implement the Yule-Walker equation by a combinational 
circuit, an algorithm that has no conditional branching must 
25 be found, and this is an essential object for the algorithm 
of the invention. 

In this case, the Peterson method known for decoding of the 
Reed-Solomon codes can be used as the algorithm of the 
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invention. Using the Peterson method, the Yule-Walker 
equation can be solved directly f and the solution of the 
Yule-Walker equation can be represented as determinants by 
the Cramer formula: 



10 



Y(J) 
A 0 = 



A 2 = 



~(i) 

^~ i=l 



So 
Si 



Si 

s 2 



Si-i Si 
So Si 



Si-i 
Si 

S21-2 



S1-1-1 Si-j 

Si-i+i Si-i+2 ' 

Si Si+i 



S1-1 

S21-2-2 
S21-1 

S21-1 



Si s 2 
s 2 s 3 

Si Si+i 



Si 
Sj+i 

S21-1 



= 1, .... 1-1 



Therefore, determinants Ao } need be obtained for each 1=1, 

Me) 

t, and the determinants A A , i = . . . , e need be 
calculated for the number of errors e. 



However, when the expansion of the determinants is provided 
by using a circuit, the required number of multipliers is 
dramatically increased as t is increased, so that it is 
difficult for the determinants to be directly expanded. 
15 Therefore, in this invention, the number of calculations is 
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reduced by the recursive structure of the Hankel matrix. 
The calculation of Af t(J) using the Katayama-Morioka method 
will be explained. 

5 When the calculation algorithm for A ( p in the 
Katayama-Morioka method is written for 1 = 1 to 1 = 4, the 
form in Fig. 14 is obtained. 

For comparison, a method devised by Koga will now be 
10 described as another method for recursively calculating the 
Hankel matrix. According to the method by Koga, new error 
locator polynomial 1>1+ 2uD(X f Y) is defined, 
wherein 



;1+2u D(X, Y) = 



( Si + YX 1 S i+1 + YX i+1 ••• S i+u + YX i+ » ^ 
S i+1 + YX i+1 S i+2 + YX i+2 S i+u+1 + YX 1+ » +1 



^S i+u +ra i+u S 1+v+l + YX 1+u+l ••• S 1+2u +YX^ 



J 



15 

To calculate this new error locator polynomial, the Hankel 
matrix i;i+2 uQ, which has the i-th syndrome S x as the (1, 1) 
element, is employed, and a determinant obtained by 
symmetrically removing multiple rows and columns from this 

20 Hankel matrix is defined as a Q determinant. When the 
subscript numbers of syndromes that appear as a diagonal 
element are designated in order beginning at the upper left, 
only one Q determinant can be determined. In this case, the 
Q determinant can be represented by the row of subscripts 

25 [ai, a 2 , . . ., a p ] . According to the Koga method, an 
algorithm is presented by to calculate error locator 
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polynomial 1;i+2u D(X, Y) the Q determinant. 



o 



All of the above described conventional methods have the 
following problems. First, for the algorithm used in Fig. 
5 14 for calculating A$ J) , new terms sequentially appear on the 
right due to the asymmetry of determinants to be calculated, 
and as a result, as the size of a matrix is increased, the 
number of multipliers that is required is increased 
accordingly. Thus, an algorithm for which the combinatorial 
10 increase of the number of multipliers is as small as 



y possible is preferable. 



% As for the Koga algorithm, the Q determinant defined by Koga 

|p is symmetrical, and a reduction in the number of multipliers 

s ' n is carried out. However, there is a limitation on the use 

of the Koga algorithm; this algorithm can be applied for BCH 
g codes or Reed-Solomon codes only when the minimum distance 

is an even number. Although it is disclosed in the Koga 
pj algorithm that this limitation can be eased, such an 

20 application example is limited to the binary, narrow sense 

BCH code. 



According to the invention, since the combinational circuit 

is applied for an optical communication system, it is 

25 required as an object that the decoding of the (255, 239) 

Reed-Solomon code (minimum distance = 17) be efficiently 

carried out by the combinational circuit. Therefore, an 

algorithm is needed that can perform an efficient 

calculation by a specific method, regardless of whether the 
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minimum distance is an odd or even number. 

£L: Definitions of terms used for this invention 

Before a detailed explanation is given for the algorithm of 
5 this invention, definitions for the terms used for this 
invention will be given. 
(1) Syndrome 

Generally, when a primitive element of the Galois extension 
field GF(2 m ) is defined as a, and h < 2 m -l is a positive 
10 integer, a 2 m -element cyclic code that has a code length of n 
= 2 m - 1 and that employs 
G(x) = (x-1) (x-a) (x-a 2 ) . . . (x-a 11 " 1 ) 

as a generator polynomial is defined as the Reed-Solomon 
codes. That is, when k = n - h and when the k-th degree 
15 polynomial having k information symbols as its coefficients 
is defined as M(x), M(x) and x n ~ k are multiplied, and the 
result is divided by G(x), as follows, to obtain the residue 
R(x) . 

20 M(x)x n - k = Q(x)G(x) + R(x) 

Then, the polynomial (transmission polynomial) that has as 
coefficients a coded sequence having the length n is defined 
as 

W(x) = M(x)x n - k - R(x) = Q(x)G(x). 
25 At this time, the coded transmission sequence is represented 
as systematic code, of which k information symbols are 
located on the left and h = n-k check symbols follow these 
symbols. The minimum distance d min = h + 1 of the 
Reed-Solomon code, and the maximum number of correctable 
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errors t = [h/2] are provided. 

The following decoding algorithm is given for making an 
estimation of the original transmission sequence based on a 
5 received sequence. 

(2) Calcu lation of syndromes and detection of an error 
Assume that errors have occurred, and that the locations of 
the errors are denoted by i 0 , . . . , i H and the error values 
are denoted by E lQ , . . . , E 2>1 . A polynomial having E ±0 , . . . , 
10 as coefficients is defined as 
E(x) = E u x x °+. . .+E il _ 1 x 1 ^ , 

a polynomial having a received sequence of b 0 , . . . , b n -x as 
coefficients is provided as 

Y(x) = n ib i x 1 = W(x) + E(x) 
15 where Y(x) is defined as a reception polynomial. 

Then, the reception polynomial Y(x) is employed to calculate 
the following syndromes: 

Si - Y(a l ) e GF (2 m ) , i = 0 , 1 , . . . , 2 1 - 1 . 

In this case, since W(a x ) = 0, i = 0 , 1, . . . , 2t-l is 
20 established, and the obtained syndromes satisfy S L = Eta 1 ), i 
= 0/ If • • -a 2t-l. Therefore, when there is no error, all 
the syndromes are 0, so that the value of syndromes can be 
used to determine the presence of errors. 

25 121 Determination of the number of errors and the location 

of errors 

Assume that the number of errors that occur is 1 and that 
the location of the error is i 0 , . . . , i^, i.e., assume 
that the values of Jb l0 , . . . , 2^ are incorrect. In order to 
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determine the number of errors 1 and the error locations 
represented by i 0 , . . ., i W( the following polynomial, 
having a~ io , . . . , a"- 1 " as roots, is defined: 

A»(x) = n(l-a i *x)=l+A? ) x+ . . . +A<- I V- 1 +A ( i V 

5 where a" io , . . . , a"^- 1 are called error locators and A (1) (x 
called an error locator polynomial. 



is 



Further, . . . , A ( / } are coefficients of the error locator 

polynomial and are provided using the elementary symmetric 
10 function of a u , . . . , a 11 - 1 . 



The unknown quantities A^, . 
simultaneous linear equations 



, A ( / } satisfy the following 



So 
Si 



Si 

s 2 



Si-i 
Si 



\ Si-i Si 



s 2 



1-2 J 



^ A? ) 



s, 



^ 



\ 521-1 J 



15 This is nothing but the Yule-Walker equation explained in 
section A. While at this step "1" is an unknown number, it 
is known that when the number of errors that actually 
occurred is 1 ^ e g t, the Hankel matrix on the left is 
regular when 1 = e and is irregular when t £ 1 > e. 

20 Therefore, for 1 = 1, . . . , t, only the determinants for 
the Hankel matrix on the left need be calculated, and a 
maximum integer that is not 0 can be defined as the number 
of errors e. When the above equation is solved with 1 = e, 
the error locator polynomial can be obtained. 

25 

In this invention, to specify error locations, the error 
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locators, i.e., the roots of error locator polynomial A (e) (x) 
= 0, need be calculated. For this calculation, a method can 
be employed whereby a" 1 , i = 0, 1, . . . , n - 1, is 
substituted in to determine whether the result is the zero 
5 point of the error locator polynomial. This method is 
called the Chien search method. When the zero point of the 
error locator polynomial is denoted as a' 10 f . . . , a -2 - 1 , i 0 , . . 
. , i e -i provide the actual error locations. 

10 (4) Calculation of error values 

Error values can be obtained by solving the following 
Vandermonde simultaneous linear equations: 
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In this case, polynomial S(x), having a syndrome as a 
15 coefficient, is defined as 
S(x) = S 0 + Six + S 2t -ix 2t - 1 , and 
Q(x) = A (e) (x)S(x) mod x 2 ^ 1 

is defined, wherein Q(x) is called an error evaluator 
polynomial. In this case, the solution for the Vandermonde 
20 simultaneous linear equation can be obtained by calculating 

^ ~ A'Ca" 1 *) r 1 ~ 0 t • • • r e- 1 

This is called the Forney algorithm. When the error 
locations and the error values are obtained, only these need 
be subtracted from an input digital signal, so that a 
25 digital signal for which errors were corrected can be 
output . 



53 



JP920010023US1 



D_= Algorithm to solve Yule-Walker equation of this 

invention 

The object of the present inventors is to find an efficient 
algorithm that employs the combinational circuit to obtain 
the solution for the following Yule-Walker equation, defined 
over GF (2 m ) : 



f So •■■ Si- ^ 



-i 



^ Sj-i "' S21-2 J 



1 



A (i) 



^ S2J-1 J 



where S 0 , Si, . . . , S 2t -i are the elements of a given GF(2 m ), 
and are unknown amounts. 

In this invention, by the Cramer formula the solution for 
the Yule-Walker equation is represented as the determinant 
form shown in Fig. 15, and a recursive structure is employed 
to obtain an efficient method for calculating determinants. 

In order to calculate determinants in Fig. 15, the focus in 
the invention is on the following Jacobi's formula. 
<Jacobi 1 s f ormula> 

A = (aij) is defined as an order n square matrix on a 
commutative ring having a unit element of 1, and the 
cof actor (i, j) of A is defined as Aij. When A#? is defined 
as the cofactor of minor A$ , with a set of subscripts, 

P ^ {il/ • • • f ir, (ii < . . . < i r )}/ v = {jx, . . j r , (j x 
< . . . < j r ) }, 

the following equation is established. 
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^detA^A^ 



In this invention, the following equation 

Ai iyi -Ai 2j2 -A i2jl -Ai^ =(det A)A^v" 2) 

can be employed, which holds when 

ia = {ii, i 2 }, v - { ji, j 2 } . 



will now be described using 



The calculation of Af t(i) 
Jacobi 1 s formula. 

First, Af t(i+1) is represented in the following form. 



r (i+i) 
A 0 = 



So 
Si 



Si-i Sj 
: Si+i 



By careful examination of this determinant, Af t(1) is obtained 
by removing the (1 + 1 - i)th row and the 1-th column from 
A A«t(j+i) ^ and A hat(i) is obta j_ ned by removing the 1-th row and 

the 1-th column. That is, since Ao at(i) andA^ t(1) are (1 + 1, 
1 + 1) and (1 + 1, 1 + 1-i) cofactors of A ( 0 J+1) respectively, 

A - A A 

^1+1, 1+1 - A 0 , ai+i t i+i-i - Ai+i-i, 1+1 - A i , 



while ii = ji = 1 + 1 - i, i 2 = j 2 - 1 + 1 is set in Jacobi 1 s 



formula. Further, the (1+1-i, 1+1-i) cofactor of A£ 



nat(l+l) 



IS 



defined as r ( / +1) , the structure of which is shown in Fig. 16. 
And by Jacobi' s formula, the following equation is obtained: 



-(1+3)' 



T<-0 



M-0 
+ (A- ) 



2 _ 



'<i+l) T 



r (i) i = 1 
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When Jacobi ! s formula is employed, the calculation of A^ at(i) 
results in the calculation of r ( / +1) , which is the 
determinants of symmetric matrices. It should be noted, 
however, that, to obtain Af t(i) , not only the calculation of 
r ( / +1) but also the 2x1 multiplications and the calculation 
of 1 square roots are required. Since the calculations for 
the square root and the square calculations can be performed 
as linear calculations for GF(2 m ), these calculations can be 
implemented as a circuit substantially at the similar cost 
as that of an addition. Therefore, only a very small cost 
is required, compared with a multiplier that is a non-liner 
operating circuit. Therefore, the present inventors focused 
on only the multipliers, and discussed the number of them 
that would be required. Since the characteristic of the 
GF(2 m ) is always 2 and all F ( p are symmetrical, the algorithm 
proposed here always cancels terms that are generated from 
arrangements that are asymmetrical to the diagonal line in 
the process for expanding the cof actor of determinants. For 
example, when the cofactor expansion is calculated for a 3 x 
3 symmetric matrix, the following equation is obtained 



aba 
b d e 
c e f 



■■adf+ae 2 +b 2 f+bec+c 2 d+bec=adf+ae 2 +b 2 f+c 2 d 



Since the term "bee", which is generated from the 
arrangement asymmetrical to the diagonal line, always 
appears twice, this term is canceled. Thus, when the 
algorithm of the invention is used for a combinational 
circuit including multipliers, the required number of 
multipliers can be reduced. 
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The general form of the algorithm for the recursive 



O 

o 



t + 1, i = 0, 1, 



calculation of r ( P , 1 = 1, 2, . . 
t, is provided as follows. 

5 o. - i, r ( 0 2) = So, if = s 2 

1. when 1 > 2, i = 1, 

1q - -321-41 0 ^2i-4-^l 

2 . when 1 > 2, i = 1, . . . , 1-1, 

first, one auxiliary amount for describing the algorithm is 
10 defined. 



When {ii, 
det[{ii, . 

det[{ii , . . . 



i n } is defined as a set of subscripts, 
, i n }] is defined as 



in}] = 

Sin Si 2+ i n - Xl 

15 Specifically, det[{i lf 



511 Si 2 

5 1 2 S2i 2 -ii 



Si z +i n -i 1 



. , i n }] is the determinant of a 
symmetric matrix wherein the first row is S llf ...,S in and the 
(p, q) element is S ip+itril . This determinant is obtained by 
symmetrically removing several rows and columns from the 



Hankel determinant A 



And is calculated by using 



20 det[{ii, . . . , i n }] as follows. 

rS° =s 2 i- 2 r£?+ zf sl^detKo, 1, . . . , i-2}-{i-i-i, 1-1-*}] . 

In this equation, det[{0, 1, . . 1-2} - {1-1-i, 1-1-k}] 

is the determinant for a symmetric matrix that is obtained 



by symmetrically removing, from Y\ 



(i-i) 



1-1-i and 1-1-k rows 



25 and 1-1-i and 1-1-k columns. Note that, when k = 1 and i = 
1, the determinant matches expression of r ( 2 i" 2 2) , r ( ^T 2 2) . 
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3. Generally, det[{ii, . . i n }] is calculated as follows. 
det[{ii, . . . , i n }] = 

S 2in - 11 det[{i lf . . . , + S^_ ii+ijt det[{ii, . . . , i^i}-{i k }] . 

E. Application of the algorithm of the invention for the 

decoding of the Reed-Solomon code 

An explanation will now be given for the embodiment wherein 
the algorithm of the invention for solving the Yule-Walker 
equation described in D is applied for the Reed-Solomon 
codes. Generally, it is assumed that the order of the 
Yule-Walker equation (the number of unknown quantities) is 
known. However, for the decoding of the Reed-Solomon codes, 
since the order is also unknown, this must also be 
determined . 

(1) Calculation of T{ 1} 

When a sequence of syndromes, S 0 , Si, . . ., S 2t -i, is 
provided, 

rf , 1 = 1, 2, . . ., t+1, i = 0, . . ., t 

is calculated in accordance with the algorithm explained in 
D. During this calculation, 

- r (i+1) 1-1 +- 

/v o ~ 1 o ,r — i, ...,r 

is also calculated. It should be noted that for fast 
decoding of the Reed-Solomon code, the present inventors 
have taken into consideration the fact that the algorithm is 
implemented as a combinational circuit. However, use of the 
error correction algorithm of the invention is not limited 
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n 



P 



to the combinational circuit; it can be employed as an error 
correction device by a sequential circuit. 

(2) Determination of the number of errors 
5 Assume that the number of errors that actually occurred is 
represented as e. Based on the value of 
_ WJ+D 1-1 t 

e can be obtained as the maximum "1" that satisfies Ao at(i) * 
0. 

10 

§ (3) Determination of an error locator polynomial 

M2 When e < t is established as the result of the determination 

Pi 

,p of the number of errors, since Ao at(efl) =0, in accordance 

W with the algorithm of the invention, the above equation can 

s 15 be simplified as follows: 

Ai ^V r i A o , 2 = 1, . . . , e . 

Since the error locator is the zero point of the error 
locator polynomial, the error locators are values unchanged 
by multiplication of the coefficients of the error locator 
20 polynomial by a constant. Therefore, the following quantity 



can be used instead of A* at(e) . In other words, the 
multiplication appearing in the above equation is not 
required. When e = t is established, the error locator 
25 polynomial is calculated in accordance with the following 
equation : 

Ai =vn 'a 0 +a 0 ry x , 2=1, . . ., e. 

At this time, according to the algorithm proposed by the 
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present inventors, in appearance, syndrome S 2t , which can not 
be calculated, seems to be necessary when the minimum 
distance is an odd number (== 2t + 1) . However, since the 
equation of this invention serves as the identity of the 
syndromes, it can also serve as the identify of the syndrome 
S 2t . Furthermore, since Af t(t) does not include syndrome S 2t , 
when syndrome S 2t appears during the cofactor expansion for 
Af £(t+1) and T^ t+1) , it should always be canceled. Specifically, 
since a term that includes S 2t and that appears during the 
cofactor expansion of r< t+1) is T i -! 1 S 2t , when rf +1) Af ^ is 
expanded, a term including S 2t is Ao at(t) r£\ 5 2t . Further, since, 
of the terms that appear during the cofactor expansion of 
A J*t<fcfi) f Ao at(t) S 2t includes S 2t , a term including S 2t when A ( 0 t+L) r^i 
is expanded is Ao at(t) r ( i l ) 1 S 2t - Therefore, all these terms must 
always be canceled . 

It is therefore understood that, of the terms that appear 
during the cofactor expansion of A ( 0 C+1) and r ( / +1) , terms having 
S 2t as a coefficient need not be calculated. In this manner, 
the algorithm of the invention can be applied for 
Reed-Solomon codes having an arbitrary minimum distance. 
Further, since the multiplication of the term including S 2t 
is not necessary, from the viewpoint of the reduction in the 
number of multipliers, the algorithm of the invention is 
superior to the Koga algorithm. In addition, as is 
described above, the calculation of a square root can be 
implemented as a circuit at the same cost as that for 
addition, and only a small cost is required compared with 
the cost of a multiplier. 
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F. Example of application of the error correction algorithm 
of the invention for the decoding of Reed-Solomon codes 
An explanation will now be given for a case wherein the 
error correction algorithm explained in E. is employed for 
the decoding of Reed-Solomon code for t = 4. When t = 4, 
according to the invention the following equations are 
determined. It should be noted that for simplification the 
determinant is represented as det[{i x , . . . , i n }] = detii . 



(1) Calculation of IT , i = 0, . . 1-1, .... 5 

The calculation results obtained by the invention are shown 
in Fig. 17. 

(2) Determination of the number of errors 

Since = r ( 0 i+1) , 1 = 1, . . . , 4 is determined by r* 2) , 

obtained in (1) , and the number of errors e can be 
determined as the maximum 1, 1=1, 2, 3, 4, that satisfies 



(3) Determination of an error locator polynomial 
When, for example, e = 2 is ascertained using the 
calculation in (2), the error locator a io , a 11 can be 
obtained by solving the following algebraic equation: 



When e = 4 is ascertained, as is described above, the error 
locator can be obtained by 



n • 



A 0 * 0. 
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It should be noted that, as is described above, the term 
including the syndrome S 8 need not be calculated when the 

rf } , r ( 0 6) = A ( 0 5) calculations are performed. 

Fig. 18 is a schematic flowchart for the error correction 
algorithm of the invention. In the error correction 
algorithm of the invention, first, at step 200 syndromes S 0 , 
S 2t -i are input, and at step 201 an error locator 
polynomial T is calculated. When r ( 0 2) , . . . , T ( 0 t+1) are obtained, 
at step 202 the number of errors is determined to be the 
maximum integer m that satisfies Ao at(in) = * 0. Then, at 

step 203, a check is performed to determine whether the 
number of errors e is equal to the maximum number of errors, 
and when e = t (yes), at step 204 an error value is 
calculated using rf 1) =A^ t(e) / . . . , rt^, rf^Af t(efl) . Whereas when 
e ^ t (no) , at step 205 an error value is calculated by only 
rr } =Ao at(e) , . . . , r<r 1} , and at step 206, Af t(e) , . . . , A?** is 
obtained. 

G. Calculation circuit when the algorithm of the invention 
is applied for the calculation of an error locator 
polynomial 

Fig. 19 is a block diagram showing a circuit for calculating 
an error locator polynomial based on the algorithm proposed 
by the invention. Fig. 20 is a diagram showing a circuit 
for calculating an error locator polynomial by the algorithm 
of the invention. This circuit comprises a {Tf } } 

calculation block 100, a circuit block 102 for calculating 
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the number of errors and a circuit block 104 for determining 
an error locator polynomial. 

The functions of the blocks in Fig . 19 will now be 
described. A series of syndromes that a sequential circuit 
has obtained using the input digital signal are transmitted 
to the circuit block 100. In the circuit block 100, these 
syndromes yield rf , m = 1, 2, . . t+1, i = 0, . . t, 
in accordance with the algorithm of the invention. This 
corresponds to (1) for the detailed explanation of the 
algorithm. 

Following this, the circuit block 102 employs the obtained 
value , m = 1, 2, . . . , t + 1 to calculate the number of 
errors e, and outputs , i = 0, . . . , e, which 

corresponds to the value of e. When e = t, in addition to 

the above, r ( 0 t+2) =A ( 0 t+1) is also output. This corresponds to 
detailed explanation (2) for the algorithm. The circuit 
block 104 then employs r^ 1} , i = 0, . . . , e to calculate 
the coefficients of the error locator polynomial. This 
calculation is performed in accordance with the process 
corresponding to detailed explanation (3) for the algorithm. 

The algorithm of the invention has been used for the 
combinational circuit in order to perform fast decoding of 
the Reed-Solomon codes. However, the algorithm of the 
invention can also be used for a sequential circuit in order 
to reduce the circuit size. 
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H. Circuit size when the algorithm of the invention is used 
for the decoding of the Reed-Solomon codes 

An explanation will now be given for the size of a circuit 
when the algorithm of the present invention is used for the 
decoding of the Reed-Solomon codes. As is described above, 
the calculation of square roots and the calculation of 
squares can be performed by a circuit having substantially 
the same cost as an addition circuit, and compared with a 
multiplier, the cost required is very small. The present 
inventors have focused only on the multipliers, and 
discussed the number of multipliers that are required. 

Table 4 shows the number of multipliers required by the 
algorithm of the invention in a range extending from t = 1 
to t = 8. In Table 4, for comparison, the number of 
multipliers required for each of the conventional examples 1 
and 2 is also shown. 



Table 4 

Comparison of the number of multipliers 



Max imum 
number of 
correctable 
errors t 


1 


2 


3 


4 


5 


6 


7 


8 


Algorithm for 
conventional 
example 1 


0 


3 


17 


48 


117* 


255* 


548* 


1111* 


Algorithm for 
conventional 
example 2 


2 


9 


22 


49 


98 


189 


351 


640 


Algorithm of 
the invention 


2 


7 


21 


46 


94 


179 


331 


597 



Note: a value having an appended * is an estimated one 
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As is apparent from Table 4, while taking the required 
number of multipliers into account, the algorithm proposed 
in this invention is superior in all number of errors t to 
the algorithm (conventional example 2) proposed by Koga. 
Further, the use of an algorithm for the decoding of (255, 
239) Reed-Solomon code (t = 8) is especially important for 
the optical communication field; however, the Koga algorithm 
can not be so employed because the minimum distance of 
Reed-Solomon code is an odd number (= 17) . Since the 
algorithm of the invention can be used for Reed-Solomon 
codes having an arbitrary minimum distance, it can also be 
used for (255, 239) Reed-Solomon code. This is shown in 
Table 5. 

Table 5 

Comparison of the application ranges of the Koga algorithm 
and the algorithm of the invention 



Minimum 
distance of 
code 




15 


16 


17 


18 




Koga 

algorithm 




X 


O 


X 


O 




Proposed 
algorithm 




o 


o 


O 


o 





(x indicates the algorithm can not be used, and O indicates 
the algorithm can be used. While the Koga algorithm can be 
used only for code having an even minimum distance, the 
proposed algorithm can be used for code having an arbitrary 
minimum distance. (255, 239) Reed-Solomon code standardized 
by the ITU has a minimum distance of 17.) 

The calculation algorithm in conventional example 1 
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(Katayama-Morioka) can also be used for Reed-Solomon code 
having an arbitrary minimum distance. However, from the 
viewpoint of the required number of multipliers into 
account, when t is equal to or greater than 4, the algorithm 
proposed in this invention requires a smaller number of 
multipliers than does the algorithm of conventional example 
1. It has especially been found that when t = 8, the 
algorithm of the invention can reduce the number of 
multipliers by about 50%. And as for a circuit size, when t 
= 8, 10K gates are currently required for the calculation of 
the error values. For conventional example 1 about 80K 
gates seem to be required, while the employment of the 
algorithm of the invention can reduce the gates for the 
calculation of an error polynomial to about 40K gates. 

Fig. 20 is a schematic diagram showing an error correction 
device according to the invention. The error correction 
device in Fig. 20 comprises: an encoding block 110, for 
receiving and encoding a digital signal; an input block 112, 
for receiving the encoded digital signal ID and for 
calculating syndromes; a process block 114, including a 
decoding circuit; and an output block 116, for correcting an 
error using an error location and an error value that are 
output and for outputting the resultant digital signal OD. 
The encoding block 110 receives the digital signal, which is 
transmitted by interleaved wavelength division multiplexing, 
converts the signal into Reed-Solomon codes, for example, 
and transmits the encoded digital signal to the input block 
112. The input block 112 employs a sequential circuit to 
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calculate syndromes for the received digital signal, and 
transmits the syndromes to the process block 114. 

The process block 114 includes a decoding function employing 
5 the algorithm of the invention, and calculates error 
locations and error values. The error locations and the 
error values that are obtained are transmitted to the output 
block 116, the error is corrected, and the resultant digital 
signal is output. The above described error correction 

10 circuit can be provided as an error correction device 
comprising multiple hardware components, or a semiconductor 
technique may be employed to provide a semiconductor device, 
such as an ASIC, for which the individual functional blocks 
of the error correction circuit are implemented on a silicon 

15 wafer. In addition, the algorithm of the invention can be 
mounted as firmware for the error correction device, or may 
be provided as a computer-readable program that is recorded 
on a storage medium, such as a floppy disk, a hard disk, an 
optical disk or a magneto-optical disk. The program of the 

20 invention may be written in an arbitrary object-oriented 
language or a programming language such as C, and stored on 
the above mentioned storage medium. 

As is described above, according to the present invention, 
25 it is possible to provide a program and storage medium, that 
can extremely efficiently correct errors in the fast optical 
communication field, perform the error correction method, 
decoding method and error correction. 
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Description of the Symbols 
10: Input unit 
12: Processor 
5 14: Output unit 

16: Syndrome calculator 

18: Error locator polynomial calculator 
18a: Register 

20: Error value polynomial calculator 
10 22: Register 
fa 24: AND gate 

M 26: XOR gate 

if*? 

Pi 28a, 28b: imn buffer 

.E 

h 40, 42: Multiplier 

H 15 45a, 45b, 62, 64: Output 

3 

46: Multipliers 

W 47: Adder group 

p 52, 54, 56, 60a, 60b: XOR group 

Q 60: Downstream XOR group 

20 66, 68, 70, 72, 77, 78: Input 

80, 82: AND group 

84: Addition circuit 
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